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| 1.0 INTRODUCTION 


This report describes the work performed by Radian 
Corporation on Contract No. NOQ0Q024-71-C1246 from 1 July 1971 
i to 1 September 1972, This work consisted of on-site support 

| to Naval Undersea Research and Development Center (NUC) personnel 
in the area of sonar transducer analysis and design. In addition, 
| work was performed to investigate the cffects of diffraction and 
r transmission characteristics of planar baffles on passive array 


| performance. 
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2.0 TRANSDUCER ANALYSIS AND DESIGN 


This section describes the work performed in the on- 
site support to the Transducer and Array Systems Division of NUC, 
San Diego, This work involves the computer-aided transducer 
analysis and design using SEADUCER (Steady-state Evaluation and 
Analysis of TransDUCERs). This computer model is designed to 
mcdel any linear system which can be described as an interconnec- 
tion of multi-port networks; however, its function for this con- 
tract was related strictly to sonar trans..rers, The on-sixe 


support included the following tasks: 


Reformat magnetic tapes sup: .ed to NUC 

by General Dynamics/Electric Boat Division 
that contained impedances and velocities to 
be compatible with SEADUCER input format, 


Design and set up SEADUCER rums to analyze 
the properties of the 3C5 transducer ele- 


ments varying the element configuration. 


* Review documentation of Modest Improvement : 
Plan (MIP) for the TR-208 active trans- 


ducer clement. 


wie remme 


. Convert SEADUCER from UNIVAC 1108 EXEC 8 
to UNIVAC 1108 EXEC II, and execute several 


example runs. 


QQ wee te ee 
xX 


Reformatting the 9 track magnetic tapes to a format 
compatible with the SLADUCER input format was accomplished via 
subroutines written for the NUC UNIVAC 1108 facility. ‘The data 


on the tapes supplicd by GD/EB were impedance matrices and 
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velocities for different heads for the 3C5 transducer. The 
impedances matrices were reduced to determine an array of unique 
impedances to be used in the computation of various performance 
variables as a function of frequency. This process was done ,to 
reduce the computation time associated with the performance 


variables, 


Radian assisted NUC personnel in the design and set-up 
of a single frequency design run to analyze the 3C5 element 
(Run No. 80168). The analysis consisted of examination of the 
element without a can, and then adding the can and mounts to 
determine their effect on the element performance. The 3C5 ele- 
ment was also analyzed with different transducer heads that were 
designed by GD/LB. 


The documentation of the Modest Improvement Plan was 
reviewed by Radian, This documented work performed by NUC and 
assisted by Radian for improvement at the performance of the 
TR-208 active transducer element. The analysis consisted of com- 
puted critical performance variables as a function of frequency 


and displaying these via a Calcomp plotter, 


The SEADUCER computer model was implemented at the 
Radian facility which utilizes a UNIVAC 1108 operating under the 
EXEC Il operating system. This required some minor changes in 
control cards from the NUC UNLVAC 1108 EXEC 8 vexsion. Radian 
converted the magnetic tape supplied by NUC using a program 
supplied by UNIVAC. After conversion of the tape, the model 
was checked by running, the example runs that are included in 
the SEALDUCER model documentation, Two of these example runs 


are included in Appendix A, 
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3.0 DIFFRACTION INVESTIGATIONS 


This section describes the support effort that Radian 
has provided in connection with NUC investigations on passive- 
array baffle design. These efforts have been concerned with 


a tel al ale ca a 


developing a computer model to investigate the effects of dif- 
fraction and transmission characteristics of planar baffles on 


passive~array performance. vi 


The modeling efforts represent one facet of an overall i i 
exploratory development program whose objective is to provide 
design guidelines for improved baffles to enhance the low fre- 
quency, passive performance of wide aperture arrays. Due to the 
increased size and weight requirements associated with extending 
the conventional puffs baffle concept to lower frequencies, the 
goal of the development effort is to examine possible alternatives 
to design a functional haffle which is effective in shielding the 
array from internal machinery noise and simultaneously enhances 


low frequency, passive-array performance. 


In connection with this program, Radian's support efforts 
have involved development of a computer model to characterize the 
scattering and transmission properties of planar strip baffles. 

To date a working computer model has been developed for calculat- 
ing the pressure field in the vacinity of a planar baffie which 
accounts for the combined effects of transmission and diffraction 
from a half plane. A version of this computer program is opera- 
tional on the UNIVAC 1230 at NUC and is being used in conjunction 
with an experimental investigation of scattering from rectangular 


plates of various materials for a wide range of frequencies. In 


re ee re 


addition, this model is currently being veneralized to treat the ig 
case of a planar strip baffle according to the approach outlined 


in Appendix B. 
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| These efforts represent an initial attempt, as part of 
@ continuing investigation, toward the synthesis of a valid com- 
i | posite math model to perform parametric analyses of significant 
baffle parameters on performance of passive-arrays. 
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Appendix A 


Examples of SEADUCER Run Under UNIVAC 1.08 EXEU IL. 
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This appendix shows examples of the SEADUCER model run 
under UNIVAC 1108 EXEC II. The purpose of this task was to 
demonstrate that the SEADUCER model was written in standard 
FORTRAN and could operate under different operating systems, and 
also on different types of computers. The task also made avail- 
able to potential users of the model a version that had compiied 
and executed correct under the EXEC II operating system. 


The first example presented here is a single frequency 
transducer element design, The transduccr element is described 
by defining network interconnections and specifying individual 
network types. Then for a varicty of tail lengths, ceramic 
areas, and ceramic lengths, the problem is to determine the 
number of ceramic pieces required to minimize the voltage- 
velocity control impedance, Zoe! The second example is strictly 
matric calculation for a given transducer element configuretion, 
and was included since it checks many of the features of the 
model, 
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Scattering from Translucent Baffles 


The purpose of this discussion is to outline an approach 
for treating the combined effects of transmission and edge dif- 
fraction for plane wave scattering from a planar baffle. The 
scope of the discussion includes a brief description of an 
integral formulation of the problem and some of the difficulties 
inherent in an approximate approach. 


A rigorous formulation of the problem of interest dif- 
fers from that of an opaque baffle in the requirement of character- 
izing the response of the baffle material to a distribution of 
forces on its surfaces. In effect, this requirement is equivalent 
to the determination of a Green's function or "influence" function 
for the baffle material.* 


Knowing the influence function for the baffle material, 
one can write an integral relationship between the force distri- 
bution on the baffle surfaces and the velocity field on the batfle 
surface. Since the baffle is immersed in a fluid which is assumed 
inviscid, the force at a given point on the baffle is normal to 
the surface. Consequently, the normal velocity ve) at a point 
Yr on the baffle surface depends on the normal force per unit area, 


i,e., the pressure, p, according to the relation: 


Vi) = | KG@lE)pG@'yas' , (B-1) 
S 
b 


* See Sec. IV of Ref. [B-1] for en analogous approach to 
transducer dome interaction. 


* 


nt 
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. 


| where K(r[r') is the n-n component of the influence dyadic [B-1] 
for the baffle material. The integral in Equation (B-1) is taken 
| over both sides of the baffle. 


| It is important to note that the kernel, K(r]r'), de- 
pends only on the geometrical and physical properties of the 

| baffle, i.e., it is independent of the conditions existing at 
any point of the fluid in which the baffle is immersed. Also we 

note that Equation (B-1) is a generalized impedance relationship; 

| more correctly K(r/r') has the units of an admittance per unit 
area. As such, Equation (B-1l) is a mathematical formulation of 

| the transmissive properties of the baffle material. Furthermore, 


any baffle material having linear constitutive equations gives 
| rise to an impedance relationship of the form given in Equation 


(B-1)*, i.e., the form of Equation (B-1) does not depend on any t 
t assumptions regarding baffle geometry. ' 
| To complete the specification cf the scattering of a 
plane wave 
py(r) = e (B-2) 
| by the baffle, we use the Helmholtz representation of the total ; | 
y , 
| 


field: ' 
| 
p(E) = p(t) + if Lp) se Elz") - PE ele) fas’, 
b 
| (B-3) 


——— 
nee es 


* For cases of locally reacting surfaces, the kernel is pro- 
| portional to a (surface) "delta function". 


. 7 .7 
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where p(r) is the total pressure field at the point = in the fluid, 


and g is the free space Green's function for the fluid nedion.> 


The above relationships, in conjunction with the defini- 
tion 


? (B-4) 


where k is the wave number and pc is the fluid, acoustical imped- 
ance, serve to define a solution to the problem. In particular, 
Equation (B-1), and the limiting form of Equation (B-3) as the 
field point Y approaches the baffle surface, give rise to a pair 
or coupled integral equations which, by virtue of Equation (B-4), 
relate the unknown pressure and velocity fields on the baffle 
surface, p and V_, respectively, to the incident pressure field 
Py A simultaneous solution of these integral equations defines 
the terms appearing in the integral of Equation (B-3), hence 
Equation (B-3) defines the pressure field at every point in the 
fluid. 


While the above approach is conceptually applicable to 
general baffle shapes in either two or three dimensions, the 
practical aspects of implementing the formalism is limited first 
of all to baffles for which the influence function can be evaluated 
analytically and secondly by computational difficulties of obtain- 
ing an approximate solution to the coupled integral equations.* 


In view of the difficulties mentioned, it would be 
desirable to have an alternative, approximate approach to repre- 


sent the essential features of the pressure field. In order to 


* Additionally the formalism would need to be modified for 
treating "'reasonance'' wave numbers associated with the baffle. 


B-3 
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derive such an approximation it is necessary to limit the dis- 
cussion to specialized geometrical configurations. In particular, 
we limit the subsequent discussion to the case of a planar baffle. 


As a first approximation, we assume that the pressure 
and velocity fields on the baffle surfaces are not influenced 
appreciably by the presence of the edge. This assumption cor- 
responds in some sense to the Kirchhoff approximation used pre- 
viously for analysis of opaque baffles. Furthermore, the assump- 
tion enables one to write 


eik-r + R(kyerky °F on S$ 
pe") | Re 
Sp T(k)e* = 


where R(k), T(k) denote, respectively, the (amplitude) reflection 
and transmission coefficients for a layer of the baffle material 
which incorporate the physical properties of the baffle’, ky 
denotes the propagation vector of the reflected wave, S; denotes 
the illuminated baffle surface, and Ss denotes the shadow baffle 
surface. In writing the above relationship, we have implicitly 
assumed the baffle thickness is negligible relative to the acous- 
tical fluid wavelength. This assumption simplifies the discussion 
and can easily be removed. 


* See Reference B-2 for the definition of the transmission and 
reflection coefficients for a composite layer. 


B-4 
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If, in addition to Equation (B-5), we assume an 
analogous equality for the normal derivative on the baffle 
surface, then the resulting approximate pressure field takes 
the form: : 


p@) = Fa [1-1 EI) +R@uElk) 6) 


where 
od bps [y= ont] eae 2,1 
velk) = | [ele") 32 etFF - ofr" 2a ) ids" 
S 
° (B-7) 


The above relationships result from substituting 
Equations (B-2), (B-5) into Equation (B-3) and using the fact 
that the integral over S; is the negative of the integral over 
S5 because of the assumption of negligible thickness of the 
baffle, 


As in the case for opaque baffles, the approximation 
scheme described above gives rise to an approximate Green's func-~ 
tion which violates the reciprocity principle. It has been noted 
previously (B-4) that the problem of plane wave scattering corres- 
ponds to a limiting process on the Green's function in which the 
source point moves to infinity. Furthermore physical considera- 
tions have indicated that for an opaque baffle, a more accurate 
approxination for plane wave scattering can be obtained from a 
Kirchhoff approximation to the Green's function by interchanging 
field and source coordinates before taking the limit. Whether 
this holds true for the present approximation remains to be 
investigated, but the assumption that it does gives rise to an 
approximation which replaces Equation (B-6) by an integral re- 


presentation which involves reflection and transmission coefficients 


B-5 
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integrated over all angles. This form arises by using the plane 

wave integral representation for the Green's function given in 

Reference B-2. Presumably, the latter approximation accounts 

partially for the fact that the presence of the edge modifies 

the source strengths appearing in Equation (B-3) from the zero 

order approximation given by Equation (B-5). This point re- 

quires more investigation. In the final analysis it may be re- 

quired to check any approximation with the pair of coupled inte- 

gral equations for justification of the assumptions. _ ; 


The above relationships represent approximate expressions 
for the pressure field which are applicable to either plates or 
strips. In the former case g denotes the three-dimensional free 
space Green's function and the integral is taken over the shadowed 
surface of the plate.* In the latter case, the integral can be 
considered as one-dimensional with g denoting the two-dimensional 


Green's function: 
glx!) = § Ho(klF-¥'1) (B-8) 


with H, denoting the zero order Hankel function of the first kind. 


Te eainaeae dr Teadaat eines aR 


* The functional form of the integral in Equation (B-7) is 
such that its value is determined strictly by the edge line 
of the plate. 


B-6 


Radian Cor POPATION «ase SHOAL Crete BvD. + 0.0. 8OX s8 + AUSTIN, TEXAS 70746 + TELEPHONE 612 - 44.9535 


| The subsequent discussion is concerned with a specialization of 
Equations (B-6) and (B-7) for the case of a strip baffle. 


| ao 2 
| my ns ay | 
{ 
rf 7) , 
| 


ey Direction of 
Strip Incidence 


| FIG. B-1 - GEOMETRY OF PLANE WAVE SCATTERING FROM A STRIP 


| Consider the geometrical situation depicted in Figure B-1. It = 
is convenient to reference the arbitrary field point r relative 
| to both edges in terms of two systems of cylindrical cocrdinates 


as indicated in Figure B-1., 


To evaluate Equation (B-7), we make use of the fact 
| that the Line integral can be taken over any contour having its 
end points on the edges of the strip and lying on the same side 
| of the field point as the strip. In particular, the contour can 
be taken to consist of Sy: the straight line from the left edge 
to the point a small distance 6 from the field point as indicated 
E in Figure B-1; S,: the portion of a circle of radius § centered 
at the field point and contained between rays from the field point ‘ 
[ 
| 
[ 


1 ates 


passing through the edges of the strip; Sy: the straight Line 
joining the field point to the right-hand baffle edge taken from 
a distance & away from the field point to the right-hand edge. 


B-7 
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It is important to note that the contribution from S, 
does not vanish as the radius 5 shrinks to zero. The simplest 
means of evaluating this Limit is to recognize that it is equal 
to the limit for the integral taken over an entire circle multi- 
plied by the ratio of the angle subtended by S, end 27. Recog- 
nizing that the angle subtended at Yr by S, is 7 - 0) = 85, and 
that the limit over the entire circle is we ik-r according to 
the Helmholtz integral formula, we obtain 


See es “Be nts 
a 2 ike 
Us ae U nn 3 


Also by observing that on Sy the normal derivative of G is 
identically zero we find in the limit as 6 approaches zero that 
Ue 7 ‘ dt H, (tye le , 


where we have used Equation (B-3) and the fact that k has an 
angular orientation a) relative to the left-hand edge. 


A similar result is obtained for the right-hand edge 


except that the angles are referenced with the subscript 2. Com- 


bination of the results gives 


kr 


ik-r-6, + 6, -7 sinv, , -it cosy 
wen Lee 2 1 , : 
Uc|K) = e | _*— + i dt Hite 
sin}, kro -it COS Yo > 
Z dt Ho(t)e ae (B-9) 
B-8 


where 


A method of evaluating the integrals in Equation (B-9) 
has already been discussed in the appendix of Reference B-3. An 
alternative means of evaluation will be discussed subsequently, 
but first it is convenient to consider the limiting case of Equation 
(B-9) as the right-hand edge moves to infinity keeping the field 
point r fixed. According to Figure B-1 this amounts to letting 
to approach infinity and simultaneously allowing 6, to approach 
zero. Thus 

kr 


ess ik‘r 8, - 7 sind, 1 cit cos}, 
Lim U(@r|k) = e Sh dt Hj(t)e 
rj70 0 
8470 sina, 


cope : de Hirde | 


Using the integral representation of the Hankel function and 
rotating the contour in the t- plane by 90° it may be shown that 


sino - it cosa om for -m<a<T, and 
as ' dt H.(t)e = ’ 
“0 periodic outside this 
interval. 
(B-10) 
B-9 
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This result in conjuncticn with the identity a = 7 - a, eivdent 
from Figure B-1 gives: 


a are -ik'r 
Lim Ur|k) = u(kr,, ¥,)e ; (B-11) 
a 
§570 
vhera 
Vy - 7 sind "1 -it cosy 
u(kr,, Vy) 2, oe ite ee : dt H_(tie ; 
(B-12) 


the above relation being valid under the assumption that Ay of 
lesz than 1. A similar analysis shows 


a -ik'r 
Lim U@r|k) = u(kry, vo)e : (B-13) 
ry" 
8470 


The relation (B-12) is the same as Equation (3.23) derived by a 


different meth.d in Reference B-5 if we recognize that Reference B-5 
uses che complement of Ya: It corresponds to the field diffracted 


by the edge of « black half-plane the geometric shadow occurring 
for Ws 0. 


Furthermore we have according to Equations (B-9), (8-11), 


(B-12), and (B-13) that 


-ik:r 


UG fk) oar [1 + u(kr,, ¥) + u(kry, Vo)! > 


which in conjunction with Equation (B-7) shows that the field is 
the superposition of the effects from the two edges. 


B-10 
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Making use of Equation (B-10) it follows that Equation 
(B-12) can be written as 


u(kr,¥) = £(¥) - spat i dt H,(t)e nee , (B-14) 
where 
0 y¥>o0 
£(¥) = 
-1 y< 0 


A simple approximation to the integral in Equation (B-14) can be 
obtained by using the asymptotic form of the Hankel function 
“to get 


| 1 -it cosy 
\ v(kr, y) = spas Jkr dt H, (te 
| sing go -i0r/*) . ei(l-cost)t 
"a on & ao 
Or 
y gat de SGD ER a 
v(kr y¥) ~ cos ($)|1 Ske 2 GUE | (B-15) 
where 


eat [x sin y/2. (B-16) 


Equation (B-16) is valid only for 4 = 0 and we use the 


fact that the original expression is odd to determine the approxi- 


mate representation in the shadow. Note that curves of constant 
x are parabolic. 


B-11 
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While Equation (B-16) is strictly valid only for kr 
large relative to unity, in practice the approximation compares 
well with the previous representation given in Reference B-3 for 
distances greater than 2 or 3 wavelengths. Furthermore a slight 
modification of the previous analysis shows that applying the 
approximation in Equation (B-16) to the case of a black half- 
plane gives a pressure field in the shadow region that differs 
from the Sommerfeld function Boe by the factor cos ¥/2. 


The above relationships provide an approximate repre- 
sentation for the scattering from a translucent strip baffle 
which is straightforward from a computational point of view. 
Furthermore it is not difficult to show that Equation (B-6) 
reduces to the solution for an infinite planar baffle as the 
edges of the strip move to infinity. This fact can be obtained 
most simply by recognizing that for this case the integral in 
Equation (B-7) can be evaluated in terms of the Helmholtz inte- 
gral formula if one recognizes that the integral over a large 
semicircle drawn in the shadow region vanishes. Consequently 


we find 
U(r|k) = - eik-r in the shadow 
Oo. in the illuminated region 
u(r |k,) = 0 in the shadow 
ve 
eth Z in the illuminated region. 


Substitution of these relations into Equation (B-6) gives the 
usual results for an infinite planar baffle. 
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(B-3) , (B-4) concerning dif frac- 


| Previous investigations 
tion from opaque baffies have shown that the diffracted field 

| computed according to the above approach gives good agreement 

with known exact results for strip baffles greater than or equal 

{ to four wavelengths.* Consequently one might expect similar 
criteria to apply also in the present case. On the other hand, 

| it should be noted that the above approximation does not account 
for the additional structural modes of vibration resulting from 

| the boundary conditions on the plate edges which are not present 

. in an infinite plate. The extent to which these modes affect the 
| pressure field in the vicinity of a finite plate remains to be 


investigated. 

‘ 
| 
| 
| 
| 
| 
: * Note that an opaque baffle can be characterized in terms of 

t Equation (B-6) by setting the transmission coefficient to zero 


and the reflection coefficient to plus or minus unity. 


f B-13 


i 
i 
i 
i 
i 
t 
i 


REFERENCES 


B-l. R. F. Pohler, "Concerning the Influence of 


Transducer Element and Dome Models on Transmit Receive Reciprocity", 


TRACOR Document 69-915-U, 17 October 1969. 


B-2. L. M. Brekhovskikh, Waves in Layered Media, 
Academic Fress, New York, 1960. 


B-3. R. F. Pohler, "An Assessment of the Validity of 
the Kirchhoff Approximation with Respect to the Radiation Charac- 
teristics of an Array on a Strip Baffle", TRACOR Document T70-AU- 
7442-U, 1 August 1970. | 


B-4, R. F. Pohler, "The Shielding Characteristics of 
A Surface Ship Ster:: Baffle’ TRACOR Document T71-AU-9510-U, 
20 January 1971. 


B-5. B. Baker and E. T. Copson, The Mathematical Theory 
of Huygen's Principle, Oxford Press, 1950. 


B-6. A. Sommerfeld, Optics, Academic Press, 1964. 


B-14 


i 
, 
? 
i 


. ® 
Radian Corporation 9500 SHOAL CREEK BLVD. © 9. ©. BOX P48 © AUSTIN, TEXAS 28768 © TELEPHONE $12 - 456-1525 


DISTRIBUTION LIST 
Contract No. NO0N024-71-C-1246 


Item No. 0002 


A0O1 AB Final Technical Report 


Sent to: 


Commander 

Naval Ship Systems Command 
Department of the Navy 
Washington, D.C. 20360 


Attn: 


OOV2 
PMS-385 
PMS ~-386 
PMS -387 
PMS-388 
OOV3 
OOV12 


Director 


Naval 


Research Laboratory 


Washington, D.C. 20390 


Attn: 


Sam Hanish, Code 8150 
Dr. Paul Smith, Code 6431 


Director, 


Naval 


Research Laborator- 


Underwater Sound Refere ce Division 


P.O. 


Box 8337 


Orlando, Florida 3zuv6 


Commander , 


Naval. 


Ship Research & Development Center 


Washington, D.C. 20007 


Attn: 


C. Chatten 


Commander, 


Naval 


Ship Research & Development Lab. 


Panama City, Florida 


Attn: 


Rufus Cook, Code 722 


No. Copies 


bh ee pe 


er 


SS A amr 


* 


en 


‘Magia RH 


oie 


DISTRIBUTION LIST CONTINUED 
Contract No. NGOO24-71-C-1246 


Sent to: 


Commander , 


Naval Undersea Research & Development Center 


San Diego, California 92132 

Attn: John Hickman, Code 601 
David Carson, Code 601 
Dr. Gordon Martin, Code 601 


Commanding Officer 

Navy Underwater Systems Center 

New London Laboratory 

New London, Connecticut 06320 

Attn: Ed. Parssinen, Code 2220 
B. B. Burnham, Code 2330 
Sidney Bayon, Code 2113 
Richard Strang, Code 2134 


Applied Research Laboratury 
The University of Texas 

P. 0. Box 8029 

Austin, Texas 78712 


General Dynamics Corporation 

Electric Boat Division 

Eastern Point Road 

Groton, Connecticut 06340 

Attn: George Sefcik, N09024-69-C-1344 


General Electric 

Farrell Road Plant, Building 1 
Syracuse, New York 13201 

Attn: Steve Auyer, NO0024-69-C-1006 


Raytheon Conpany 

Submarine Signal Division 

P. O. lox 360 

Portsmouth, Rhode [sland 02871 

Attn: Norman Scrotta, NQQ024-69-C-1006 


Radian Corporation @500 SHOAL CREEK BLVD. © P.O. BOX FHE © AUSTIN, TEXAS 20766 © TELEPHONE Siz - 454.9535, 


rR Re 


ho pt 


DISTRIBUTION LIST CONTINUED 


Contract No, NOQ024-71-C-1246 


Sent to: 


Stanford Research Institute 
333 Ravenwood Drive 

Menlo Park, California 94021 
Attn: Vincent Salmon 


Tracor, Incorporated 
Austin, Texas 78721 
Attn: William Moyer 


Comnander, 

Naval Ship Systems Command 
Washington, D.C. 20360 
Attns 2052 


Defense Documentation Center 
Cameron Station (Bldg. 5) 
5010 Duke Street 

Alexandria, Virginia 22314 


" Radian Corporation tsa soni crsex avo. + 9. 0. 80x oa + AUSTIN, TEXAS 70746 © TELEPHONE S12 - 456.9595 


No. Copies 


TTR at RPE Me RNS I 
t 


S 1 nee oa 
po Sans i B 


_ aneillin 


wer 


